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Axion-photon conversion induced by intergalactic magnetic fields has been proposed as an expla- 
nation for the dimming of distant supernovae of type la (SNe la) without cosmic acceleration. The 
effect depends on the intergalactic electron density n e as well as the B-field strength and domain 
size. We show that for n e ^$ 10 -9 cm -3 the same mechanism would cause excessive spectral dis- 
tortion of the cosmic microwave background (CMB). This small-n e parameter region had been left 
open by the most restrictive previous constraints based on the dispersion of quasar (QSO) spectra. 
The combination of CMB and QSO limits suggests that the photon-axion conversion mechanism can 
only play a subleading role for SN la dimming. A combined analysis of all the observables affected 
by the photon-axion oscillations would be required to give a final verdict on the viability of this 
model. 
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Supernovae of type la (SNe la) at redshifts 0.3 < z < 
1.7 appear fainter than expected from the luminosity- 
redshift relation in a decelerating Universe 0. 0. 
On the other hand, the cosmic microwave background 
(CMB) anisotropy and large-scale structure observations 
suggest that the Universe is spatially flat, with a matter 
density of approximately 30% of the critical density 0,13. 
The "concordance model" thus implies that the Universe 
must be accelerating today because it is dominated by 
a "dark energy" component (about 70% of the critical 
density) with an equation of state w = pj p ks — 1. 

The lack of a satisfactory fundamental explanation for 
this component has triggered wide-ranging theoretical in- 
vestigations of more or less exotic scenarios [6|. Some 
years ago Csaki, Kaloper and Terning (CKT I) sug- 
gested that the observed achromatic dimming of distant 
SNe la may be a consequence of the mixing of photons 
with very light and weakly coupled axion-like particles 
in the intergalactic magnetic fields. Though still requir- 
ing some non-standard fluid (e.g. with pj p ~ — 1/3) to 
fit the flatness of the universe, this model seemed ca- 
pable to explain the SN dimming through a completely 
different mechanism without apparently affecting other 
cosmological observations. 

Later it was recognized that the conclusions of CKT I 
can be significantly modified when the effects of the in- 
tergalactic plasma on the photon-axion oscillations are 
taken into account Assuming an electron density 
n e w ribaryons = Ti-yTj ~ 10 _7 cm -3 , the model is ruled out 
in most of the parameter space because of either an exces- 
sive photon conversion or a chromaticity of the dimming. 
Only fine-tuned parameters for the statistical properties 
of the extragalactic magnetic fields would still allow this 
explanation. On the other hand, Csaki, Kaloper and 
Terning (CKT II) criticized the assumed value of n e 
as being far too large for most of the intergalactic space, 



invoking observational hints for a value at least one or- 
der of magnitude smaller. For n e < 2.5 x 10~ 8 cm~ 3 , 
the photon-axion mixing hypothesis works even better 
with the plasma, because the constraints from CMB 
anisotropics via photon-axion conversion can be relaxed. 

If distant SNe la are dimmed by this mechanism, the 
same would apply to other sources. In particular, one 
would expect a dispersion in the observed quasar (QSO) 
spectra. An analysis based on the first data release of 
the Sloan Digital Sk y S urvey excludes a large part of 
the parameter space [l(|, suggesting that only for n e < 
10~ 10 cm" 3 the axion mechanism is still able to explain a 
dimming by 0.1 magnitudes or more. If the QSO spectra 
had an intrinsic dispersion at the 5% level would rule out 
axion dimming exceeding ~ 0.05 mag. Future data will 
be sensitive to yet larger regions in parameter space, yet 
QSOs will never be sensitive to very low n e . 

A similar bound has been obtained by a possible vi- 
olation of the reciprocity relation between the luminos- 
ity distance and the angular-diameter distance [Til IT^ . 
However, this constraint is less robust than the QSO one 
because it is affected by poss ibly large systematic errors 
that are difficult to quantify [l3j- 

The purpose of our paper is to further constrain the 
photon-axion conversion model by studying its effect on 
the CMB spectral shape. We will show that the low-n e 
region of parameters left open by the QSO limit is ruled 
out by our new limit, leaving little if any room for the 
axion hypothesis to mimic cosmic acceleration. 

In Sec. [H] we discuss the formalism of photon-axion 
conversion and in Sec. IIHI we summarize its effect on 
SN la dimming. In Sec. IIVI we describe the constraints 
coming from spectral CMB distortions and in Sec. IVl we 
combine our new limits with those from QSO dispersion. 
Finally, in Sec. I VII we draw our conclusions and comment 
on the viability of the photon-axion conversion mecha- 
nism. 
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II. PHOTON-AXION CONVERSION 

Axions and photons oscillate into each other in an ex- 
ternal magnetic field [i^lTMITrilri| due to the interaction 
term 



.CM 



(1) 



£ a -f = — j ga-r^fivF^a — g^E • B a . 



where F^ v is the electromagnetic field tensor, = 
\^^vpaF pa is its dual, a is the axion field, and g aj is 
the axion-photon coupling (with dimension of inverse en- 
ergy) . We always use natural units with Ti = c = fee = 1 • 
For very relativistic axions, the equations of motion in 
the presence of an external magnetic field B reduce to 
the linearized form Ha 



cx . These terms are of little importance 
for the following arguments and will thus be neglected. 
The Faraday rotation term Ar, which depends on the 
energy and the longitudinal component B z , couples the 
modes A\\ and Aj_. While Faraday rotation is important 
when analyzing polarized sources of photons, it plays no 
role for the problem at hand. 

With this simplification the Aj_ component decouples, 
and the propagation equations reduce to a 2-dimensional 
mixing problem with a purely transverse field B = Br 



id z + A4 2 



a 



0, 



(6) 



with a 2-dimensional mixing matrix 



(uj - id z + M) 




0, 



(2) 



where z is the direction of propagation, A x and A y corre- 
spond to the two linear polarization states of the photon 
field, and us is the photon or axion energy. The mixing 
matrix is 



M 2 = 



A p i A a7 
A a7 A a 



(7) 



The solution follows from diagonalization through the 
rotation angle 



^iarctanf-^ 
2 VApi 



(8) 
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M = | Ay X Ayy 

29a^yB x ^Qa^y By A a 
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In analogy to the neutrino case [18J , the probability for a 
(3) photon emitted in the state A\\ to convert into an axion 
after traveling a distance s is 



where A a = — to^/2oj. The component of B parallel 
to the direction of motion does not induce photon-axion 
mixing. The quantities with i, j — x,y mix the pho- 
ton polarization states. They are energy dependent and 
are determined both by the properties of the medium 
and the QED vacuum polarization effect. We ignore the 
latter, being sub-dominant for the problem at hand 0. 

For a homogeneous magnetic field we may choose a 
coordinate system aligned with the field direction. The 
linear photon polarization state parallel to the transverse 
field direction By is denoted as Am and the orthogonal 
one as A±. Equation becomes then 



A i 



(u-id, + M) A { \ =0, 



(4) 



with mixing matrix 



M = 




(5) 



P ( 7 ^a) 



|(A||(0)|a( S ))^ 
sin 2 (2?9)sin 2 (A osc s/2) 
sin 2 (A osc s/2) 



(A osc s/2) 2 



where the oscillation wavenumber is given by 



A 2 sc = (A pl -A a ) 2 + 4A 2 



7 ' 



(9) 



(10) 



The conversion probability is energy-independent when 
2|A a7 | ^> |A p i — A a | or whenever the oscillatory term in 
Eq. 10 is small, i.e. A osc s/2 <C 1, implying the limiting 
behavior Pq = (A Q7 s) 2 . 

The propagation over many random B-field domains 
is a truly 3-dimensional problem, because different pho- 
ton polarization states play the role of A\\ and A± in 
different domains. This is enough to guarantee that the 
conversion probability over many domains is an incoher- 
ent average over magnetic field configurations and photon 
polarization states. The probability after travelling over 
a distance r>s, where s is the domain size, is 



Here, Aj 



A CM 



•II 



/2us, A a7 = g aj \B T \/2, and 



f Ajp 1 , A pl = 
4na n e /m e de- 



fines the plasma frequency, m e being the electron mass 
and a the fine-structure constant. The A™ terms de- 
scribe the Cotton-Mouton effect, i.e. the birefringence of 
fluids in the presence of a transverse magnetic field where 



1 — exp 



2s 



(11) 



with P given by Eq. @ . As expected one finds that for 
r/s — > oo the conversion probability saturates, so that 
on average one third of all photons converts to axions. 
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III. PHOTON-AXION CONVERSION AND 
SUPERNOVA DIMMING 

To explore the effect of photon-axion conversion on 
SN dimming we recast the relevant physical quantities in 
terms of natural parameter values. The energy of optical 
photons is a few eV. The strength of widespread, all- 
pervading B-fields in the intergalactic medium must be 
less than a few 10 -9 G over coherence lengths s crudely at 
the Mpc scale, according to the constraint coming from 
the Faraday effect of distant radio sources [20J. Along 
a given line of sight, the number of such domains in 
our Hubble radius is about N ~ Hq 1 /s « 4 x 10 3 for 
s ~ 1 Mpc. The mean diffuse intergalactic plasma den- 
sity is bounded by n e < 2.7 x 1CP 7 cm" 3 , corresponding 
to the recent WMAP measurement of the baryon den- 
sity 0]. Recent results from the CAST experiment [2l| 
give a direct experimental bound on the axion-photon 
coupling of g ai < 1.16 x 10~ 10 GeV -1 , comparable to 
the long-standing globular-cluster limit [ljj. For ultra- 
light axions a stringent limit from the absence of 7-rays 
from SN 1987A gives g ai < 1 x 10" 11 GeV" 1 ^ or 
even g aj < 3 x 1(T 12 GeV" 1 0. Therefore, suitable 
numerical values of the mixing parameters are 



A a7 



Mpc" 1 

Aq 

Mpc" 1 

A p i 
Mpc" 1 



0.15 g w B nG , 

VI eVJ VleV/ 

/ (1) 

-11.1 



leV 



10- 7 cm- 3 



(12) 



and 



where we have introduced 510 = <? a7 /10" 10 GeV" 
-B n G is the magnetic field strength in nano-Gauss. 

The mixing angle defined in Eq. JHJ is too small to 
yield a significant conversion effect for the allowed range 
of axion masses because |A a | 3> |A a7 |,|A p i|. There- 
fore, to ensure a sufficiently large mixing angle one has 
to require nearly massless pseudo-scalars, sometimes re- 
ferred to as "arions" 0> E3- Henceforth we will con- 
sider the pseudoscalars to be effectively massless, so that 
our remaining independent parameters are gioB n Q and 
n e . Note that m a only enters the equations via the 
term to 2 — w 2 j, so that for tiny but non- vanishing val- 
ues of m a , the electron density should be interpreted as 
n e ,eS = \n e - m 2 a m e /(Aira)\. 

The distance relevant for SN la dimming is the lumi- 
nosity distance 4 at redshift z, defined by 



dl(z) = 



C 



(13) 



where £ is the absolute luminosity of the source and T is 
the energy flux arriving at Earth 

0,0. Usually the data 
are expressed in terms of magnitudes 



m = M + 5 log 



10 



Mpc 



25, 



(14) 



where M is the absolute magnitude, equal to the value 
that m would have at — 10 pc. After a distance 
r, photon-axion conversion has reduced the number of 
photons emitted by the source and thus the flux T to 
the fraction P 7 _^ 7 = 1 — P~ f ^ a - Therefore, the luminosity 
distance becomes 



d L /{P 1 



a/2 



(15) 



and the brightness 



m ~ -log 10 (P 7 ^ 7 ) 



(16) 



Distant SNe la would eventually saturate (-P 7 ^ 7 = 2/3), 
and hence they would appear (3/2) 1//2 times farther away 
than they really are. This corresponds to a maximum 
dimming of approximately 0.4 mag. 

In Fig. n we show qualitatively the regions of n e and 
510-BnG relevant for SN dimming at cosmological dis- 
tances. To this end we show iso-dimming contours ob- 
tained from Eq. l|16fl for a photon energy 4.0 eV and 
a magnetic domain size s = 1 Mpc. For simplicity we 
neglect the redshift evolution of the intergalactic mag- 
netic field B, domain size s, plasma density n e and pho- 
ton frequency u. Our iso-dimming curves are intended 
to illustrate the regions where the photon-axion conver- 
sion could be relevant. In reality, the dimming should 
be a more complicated function since the intergalactic 
medium is expected to be very irregular: there could be 
voids of low n e density, but there will also be high den- 
sity clumps, sheets and filaments and these will typically 
have higher B fields as well. However, the simplifications 
used in this work are consistent with the ones adopted in 
CKT II model and do not alter our main results. 

The iso-dimming contours are horizontal in the low- 
n e and low-<7io-B n G region. They are horizontal for any 
.910 B n Q when n e is sufficiently low. From the discus- 
sion in Sec. ^ we know that the single-domain proba- 
bility Pq of Eq. @ is indeed energy independent when 
A osc s| < 1, i.e. for |A p i|s/2 < 1 and |A a7 |s < 1. When 
n e < few 10~ 8 cm -3 and gwB n Q < 4, we do not expect 
an oscillatory behavior of the probability. This feature 
is nicely reproduced in our iso-dimming contours. From 
Fig. ^ we also deduce that a significant amount of dim- 
ming is possible only for 510-BnG <^ 4 x 10~ 2 . 

In CKT I, where the effect of n e was neglected, a value 
m a ~ 10 -16 eV was used. In terms of our variables, 
this corresponds to n e , e S ~ 6 x 10~ 12 cm" 3 . As noted 
in CKT II, when plasma effects are taken into account, 
any value n e < 2.5 x 10~ 8 cm" 3 guarantees the required 
achromaticity of the dimming below the 3% level between 
the B and V bands. The choice B n Q of a few and 310 ~ 
0.1 in CKT I and II falls in the region where the observed 
SN dimming could be explained while being marginally 
compatible with the bounds on B and giQ. 
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FIG. 1: Iso-dimming curves for an attenuation of 0.01, 0.1, 
and 0.4 magnitudes. The photon energy of 4.0 eV is repre- 
sentative of the B-band. The size of a magnetic domain is 
s = 1 Mpc. 



IV. CMB CONSTRAINTS 

If 7 — y a conversion over cosmological distances is re- 
sponsible for the SN la dimming, the same phenomenon 
should also leave an imprint in the CMB. We note 
that a similar argument was pre viously considered for 
photon— > graviton conversion |26|. Qualitatively, in the 
energy-dependent region of P 7 ^ a one expects a rather 
small effect due to the low energy of CMB photons (u ~ 
10~ 4 cV). However, when accounting for the incoherent 
integration over many domains crossed by the photon, 
appreciable spectral distortions may arise in view of the 
accuracy of the CMB data (at the level of one part in 10 4 - 
10 5 ). For the same reason, in the energy-independent 
region, at much lower values of n e than for the SNe la, 
the constraints on gioB n c are expected to be quite se- 
vere. The depletion of CMB photons in the patchy mag- 
netic sky and its effect on the CMB anisotropy pattern 
have been previously considered in |7j. However, more 
stringent limits come from the distortion of the overall 
blackbody spectrum. 

To this end we use the COBE/FIRAS data for the 
experimentally measured spectrum, corrected for fore- 
grounds H3 ■ Note that the new calibration of FIRAS 
is within the old errors and would not change any of our 
conclusions. The N ~ 43 data points <5>° xp at different 
energies u>i are obtained by summing the best-fit black- 
body spectrum to the residuals reported in Ref . [23| ■ The 
experimental errors cr, and the correlation indices pij be- 
tween different energies are also available. In the presence 
of photon-axion conversion, the original intensity of the 
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FIG. 2: Exclusion plot for axion-photon conversion based on 
the COBE/FIRAS CMB spectral data. The region above the 
solid curve is excluded at 95% C.L. whereas the one above 
the dotted curve is excluded at 99% C.L. The size of each 
magnetic domain is fixed at s = 1 Mpc. We also reproduce 
the iso-dimming contours from Fig. 

"theoretical blackbody" at temperature T 

^VTH^fexp^/T)-!]- 1 (17) 

would convert to a deformed spectrum that is given by 
$(w,T) = $°(w,T)P 7 ^ 7 (w). We then build the reduced 
chi-squared function 

1 N 

xl(T, A) = £ AM^ 1 ^ , (18) 

where 

A$, = $f p - $°(w i) T)P y _ n (a; i ,A) (19) 
is the i-th residual, and 

°ij = PijWj ( 20 ) 

is the covariance matrix. We minimize this function with 
respect to T |30j | for each point in the parameter space 
A = (n e , 5ioSiig) 5 i- e - T is an empirical parameter de- 
termined by the xt minimization for each A rather than 
being fixed at the standard value T = 2.725 ± 0.002 K 

a. 

In Fig. [5] we show our exclusion contour in the plane of 
n e and gioB n Q. The region above the continuous curve 
is the excluded region at 95% C.L., i.e. in this region 
the chance probability to get larger values of \t is lower 
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than 5%. We also show the corresponding 99% C.L. con- 
tour which is very close to the 95% contour so that an- 
other regression method and/or exclusion criterion would 
not change the results very much. Within a factor of a 
few, the same contours also hold if one varies the domain 
size s within a factor 10. 

Comparing our exclusion plot with the iso-dimming 
curves of Fig.^we conclude that the entire region n e < 
1CP 9 cm~ 3 is excluded for SN dimming. 

A few comments are in order. Intergalactic magnetic 
fields probably are a relatively recent phenomenon in the 
cosmic history, arising only at redshifts of a few. As a 
first approximation we have then considered the photon- 
axion conversion as happening on present (z = 0) CMB 
photons. Since P 7 _+ 7 is an increasing function of the pho- 
ton energy w, our approach leads to conservative limits. 
Moreover, we assumed no correlation between n e and the 
intergalactic magnetic field strength. It is however phys- 
ically expected that the fields are positively correlated 
with the plasma density so that relatively high values of 
ffio-SnG should be more likely when n e is larger. Our 
constraints in the region of n e > 10~ 10 cm -3 are thus 
probably tighter than what naively appears. 

V. QSO CONSTRAINTS 

Our limits are nicely complementary to the ones ob- 
tained from the effects of photon-axion conversion on 
quasar colors and spectra [10j. In Fig. [21 we superimpose 
our CMB exclusion contours with the schematic region 
excluded by quasars |3l|. The region to the right of the 
dot-dashed line is excluded by requiring achromaticity of 
SN la dimming 9J. The region inside the dashed lines 
is excluded by the dispersion in QSO spectra. Moreover, 
assuming an intrinsic dispersion of 5% in these spectra, 
the excluded region could be enlarged up to the dotted 
lines. Our CMB argument excludes the region above the 
solid curve at 95% C.L. 

A cautionary remark is in order when combining the 
two constraints. As we have discussed in the previous 
section, our CMB limits on photon-axion conversion are 
model independent. Conversely, the limits placed by the 
QSO spectra are possibly subjected to loop holes, since 
they are based on a full correlation between the inter- 
galactic electron density and the magnetic field strength, 
which is reasonable but not well established observation- 
ally. 



VI. CONCLUSIONS 

We have examined the conversion of CMB photons into 
very low-mass axions in the presence of intergalactic mag- 
netic fields. The resulting CMB spectral deformation ex- 
cludes a previously allowed parameter region correspond- 
ing to very low densities of the intergalactic medium. Our 
new limits are complementary to the ones derived from 
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FIG. 3: Exclusion plot for photon-axion conversion. The re- 
gion to the right of the dot-dashed line is excluded by re- 
quiring achromaticity of SN la dimming. The region inside 
the dashed lines is excluded by the dispersion in QSO spec- 
tra. Assuming an intrinsic dispersion of 5% in QSO spectra, 
the excluded region could be extended up to the dotted curve. 
Our CMB argument excludes the entire region above the con- 
tinuous curve at 95% C.L. 



QSO dispersion which place serious constraints on the 
axion-photon conversion mechanism. As a result, it ap- 
pears that this mechanism can hardly play a leading role 
for the apparent SN la dimming. 

The axion-photon conversion hypothesis has also been 
advocated to explain trans-GZK cutoff events in Ultra 
High Energy Cosmic Rays (UHECRs) [2^. In principle, 
UHECR photons, produced in cosmological sources far 
away, could drastically reduce energy losses while prop- 
agating in the intergalactic medium as axions. Some of 
these particles would eventually convert back to photons 
within a few GZK radii, thus justifying the observations 
of extremely high energy events as well as their isotropy. 
While one can not rule out the possibility that some UHE 
"photon-like" events at energies £>4x 10 19 eV might 
be due to this mechanism, our bounds imply that it can 
play only a subdominant role. Moreover, photons any- 
way are disfavored as candidates for the majority of the 
UHECRs. 

In summary, the CMB constraints together with pre- 
vious limits suggest that the fascinating mechanism of 
photon-axion conversion in the intergalactic magnetic 
fields does not play an important role for either the phe- 
nomenon of SN la dimming or for UHECR propagation. 
A definitive verdict would probably require a common 
analysis of SN la dimming, QSO spectra, and the Fara- 
day effect of distant radio sources, based on mutually con- 
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sistent assumptions about the intergalactic matter den- 
sity and its distribution, the intergalactic B-ficld strength 
and its distribution and correlation with the electron den- 
sity, and the redshift evolution of these quantities. Our 
results show that the low-n e escape route from the QSO 
limits is definitely closed. 
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